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Abstract: We obtain exact analytic solutions for a typical autonomous dynamical system, related to 
the problem of a vector field nonminimally coupled to gravity. 


Gravity is presently the only interaction that is not inserted in any of the unification schemes. Due to 
this fact, many toy models have appeared coupling gravity to fields that could have played significant role in 
the Early Universe, in such wise as to display desirable properties. Special interest has been put in gravity 
coupled nonminimally to other fields which could generate, among a number of new effects, a non-singular 
universe. 

In this context, dynamical systems techniques have been applied to solve the problem of coupling of 
gravity to a vector field, whose Lagrangian can be written as ^ 

/: = (f - , (1) 

where g = deig^^, {g^i, is the metric tensor), R is the scalar of curvature, G = k/Sir is the Newtonian 
gravitational constant in units that h = c = 1,\ (3 \= is an arbitrary vector field and F^i, = 

where the comma represents ordinary differentiation and the square brackets represent the skew-symmetric 
part of 

The Lagrangian (1) leads to a set of equations of motion that can be transformed into an autonomous 
dynamical system ^. By choosing the Robert son-Walker metric 

ds^ = dt^ — S‘^{t) [dy^ + + sin^9d(j)^)] , 

with the ansatz W'^ = lU^(t), the Lagrangian (1) leads to equations of motion which has a solution given by 
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'‘"w=K' -1)' 

S{t) = (<2+g2)l/2^ 

where k is the Einstein constant and Q is also a contant. Afterwards one sets X = 3{S/S) and Y = (fi/n), 
where ft = (1/fc) + (/3VE^), to write the equations of motion as 


i: = --x'^ + XY, 
3 


Y = _ XY. 


( 2 ) 


Hence, the correct interpretation of functions X and Y is important to the knowledge of the evolution 
of the model. 

Up to now there was no exact solution for this dynamical system and the analysis was carried out by 
using qualitative procedure. Here we exhibit the explicit solution in a very simple way. 

One can rewrite eqs. (2) in polar coordinates, p = (AT^ + and 0 = arctan ^, as 


p = cos^ 9 + sin^ 9 + cos 9 sin^ 9 — cos^ 9 sin , 

(3) 

• f COS 0 \ 

9 = —2/9COS0sin0 ( —^- hsin^j. 

The phase diagram related to this autonomous dynamical system can be obtained through 

9 d9' 

which leads to 


p = Cexp Q. 

where C is an arbitrary positive constant of integration and 

i (cos^ 0 + 3 sin^ 9) + cos 9 sin^ 9 — cos^ 9 sin 9 


I = 


i cos2 9 sin 9 + cos 9 sin^ 9 


d9. 


(4) 


(5) 


After a little algebra, involving only basic trigonometric relations, one can show that the solution to 
integral I readily leads to 
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p = c. 


|tan6*| 


1/2 


|cos0 + 3sin0| 

We point out that the precise knowledge of the function p{9) leads to the sketch of the desired phase dia¬ 
gram associate with the autonomous dynamical system without any qualitative analysis about the dynamical 
system in regions around the origin 

On the other hand one can choose to study the dynamical system qualitatively by means of the projection 
on the Poincare sphere This will provide usefull informations concerning the asymptotic behaviour of the 
system. As defined in Ref. 3, the Poincare sphere with unitary radius is placed over the xy-plane, this plane 
being tangent to the south pole of the sphere. Another frame is in order and it is placed in the centre of the 
sphere. Projections onto the sphere are taken by joining points of the diagram to the centre of the sphere. 
This process will gives rise to a drawing on the sphere which is projected orthogonally on the xy-plane. The 
final portrait of the phase diagram is, then, in a circle with unitary radius where the behaviour at infinity is 
identified with the border of the circle. 

Applying this method ^ to (2) one can obtain the topologies around the singular points at infinity, as in 
Fig.l. The arrows show the evolution in time. Note that besides the equilibrium points Z?(0,0), £>'(0,0), 
(7(0,0) and (7'(0,0), there are the points i3(—1/3, 0) and i3'(l/3,0), the primes indicating antipodes points. 
The system refuses to give informations on the topology around the origin by the method of linearization. 
Dulac’s test shows there are no limiting cycles. Note that the above mentioned equilibrium points are 
consistent with eq. (6). 

The complete phase diagram which compactifies infinities is shown in Fig.2. Its shape is very close to 
that in Ref.4; the difference is due to the unusual dimension used. The phase diagram in Ref.l is obtained 
in a different way: the projection is taken on the plane which corresponds to A = 0 and crosses the poles of 
the sphere, the xy-plane being tangent to the north pole of the sphere. 

Notice that any autonomous planar dynamical system 


A = P{X,Y), 

Y = Q{X,Y). 

allows a similar change of variables, and in polar form the phase diagram arises from 

dp P{pcos9, psinO) + tan9Q{pcos 9, psin9) 
d9 ^Q{pcos9, psin9) — tan9P{pcos9, psin9) 

For this reason, exact solutions may be achieved if it is possible to rewrite (7) as 

I = 
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Recently some cosmological models has been investigated in the so-called multidimensional scenario 
and in a particular case a Maxwell held in higher dimensions is coupled to Einstein-Hilbert Lagrangian. 
This also leads to a qualitative dynamical system analysis, and, as in eq. (7), the search for exact solutions 
related to this problem is in progress. 

Thanks are due to Professor C. A. P. Galvao for reading the manuscript and for helpful suggestions. 
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Figures Captions 


Figure 1: Equilibrium points at inhnity and the topology imposed by them. Empty balls means unstable 
points, and full balls stable ones. 

Figure 2: Final aspect of the phase diagram showing the equilibrium points at inhnity. Only one integral 
curve is chosen to each region in the diagram to avoid it to become entangled. 
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